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Abstract 

The inequivalence of thermodynamical ensembles related by a Legendre 
transformation is manifest in self-gravitating systems and in black hole ther- 
modynamics. Using the Poincare's method of the linear series, we describe 
the mathematical reasons which lead to this inequivalence which in turn in- 
duces a hierarchy of ensembles: the most stable ensemble describes the most 
isolated system. Moreover, we prove that one can obtain the degree of sta- 
bility of all equilibrium configurations in any ensemble related by Legendre 
transformations to the most stable if one knows the degree of stability in the 
most stable ensemble. 
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1. Introduction 



In self-gravitating systems' 1-5 ! or in systems with highly degenerate spec- 
tra such as strings' 6 ' or black holes' 7-14 ', it is well known that thermodynam- 
ical ensembles are generally inequivalent. This means that two ensembles 
related by a single Legendre transformation (or a Laplace transformation in 
Statistical Mechanics) become unstable in different situations. For instance 
there exist stable situations in microcanonical ensembles with negative spe- 
cific heat. Surprisingly there also exist situations in which canonical ensem- 
bles with positive specific heat are unstable. 

Thus, evidently, the heat capacity alone does not control the stability 
of those ensembles and one needs more reliable criteria. These can be pro- 
vided, under specific conditions (see refs. 1, 2), by the Poincare's linear series 
method. This method offers the great advantage to work without having 
to solve the eigenvalue equation which control the quadratic fluctuations 
(i.e. without having to calculate the Poincare's coefficients of stability). The 
method provides also, as we will show, a useful guide to obtain the stability 
limits of any Legendre transformed ensemble when one knows the stability 
range of one ensemble and this, whether or not these ensembles are equiva- 
lent. 

In this article we will study the origin as well as the generic properties 
of this inequivalence. Addressing the determination of stability criteria, we 
will first find the mathematical reason which leads to the breaking of the ap- 
parent "symmetry" between an ensemble and its Legendre transform. This 
is strange indeed in that if we apply a second Legendre transformation one 
is led back to the initial ensemble. This mathematical reason will also per- 
mit us to make contact with physics because the most stable ensemble will 
systematically describe the most isolated system. By isolated we mean an 
ensemble whose variable (called in this article control parameter) which char- 
acterizes it, is a conserved quantity of the hamiltonian of the system. On 
the contrary, non-isolated systems are characterized by control parameters 
kept fixed by the intervention of an external reservoir. Familiar examples 
of couples of control parameters are respectively: energy-temperature, parti- 
cle number-chemical potential, volume-pressure, angular momentum-angular 
velocity. 
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We then demonstrate the following points: 

1. The inequivalence comes from the fact that the less stable ensemble has 
a richer spectrum of fluctuations which contains one extra fluctuating 
quantity (or one more Poincare stability coefficient). 

2. From this fact one immediately concludes that when the corresponding 
eigenvalue of those extra fluctuations is positive it simply means that 
the number of unstable modes in the two ensembles is identical (This 
number can be zero in which case both ensembles are stable). 

3. More remarkable is the fact that it is always this new eigenvalue which 
encodes the changes of stability of the less stable ensemble. This means 
that whenever the more stable system approaches instability (i.e. when 
its lowest eigenvalue tends to zero) an algebraic identity dictates that 
the new eigenvalue will cross zero before the lowest one of the more 
stable system. 

Our analysis results in a qualification of the traditional criterion wherein 
the positivity of the specific heat was considered as a sufficient condition for 
the stability of the canonical ensemble. It provides also a timesaving pro- 
cedure when, being aware that ensembles might be inequivalent, one wishes 
to determine stability ranges under various conditions' 4 ' 5 ' 14 ' 15 '. As an illus- 
tration of its power, we apply our analysis to the situation (recently studied 
by Kaburaki et a/.' 10,11 ') in which the angular momentum of rotating black 
holes has a stabilizing effect. We derive results on stability that are stronger 
than those previously published. 

Section 2 

Our goal is to compare the stability limits of two ensembles related by a 
single Legendre transformation. The demonstration of the three points men- 
tioned in the Introduction proceeds as follows. We start by the analysis of 
a simplified mathematical abstraction in which all the members of the first 
ensemble are states parametrized by a single variable, denoted by x, at a 
fixed control parameter s. Subsequently the applicability of the mathemat- 
ical exercise to physical situations including many degrees of freedom (as in 
thermodynamics) is discussed. 
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In the first ensemble, equilibrium situations and stability conditions are 
provided by a potential function F(x; s)Q The Legendre transformation to 
a new potential G(x; t) with control parameter t, is given by: 

G(x;t) = F(x;S(x;t)) -tS(x;t) (1) 
where S(x; t) is the solution (unique at fixed x) of: 

t = d s F(x; s) (2) 

d designates partial derivative with all other quantities kept fixed. The equi- 
librium configurations, Xi(s) i = 1, • • -n, are the solutions of: 

d x F(x; s) = (3) 

Those linear series of equilibrium configurations coincide with the solutions 
of 

d x G(x;t) = Q (4) 

by virtue of eq(2). At this point, the equivalence of the two ensembles is 
manifest since F is the Legendre transform potential of G, with the role of 
s and t interchanged, and since they determine the same equilibrium config- 
urations. (We shall adopt the notation F-ensemble when the description of 
the states is given by F(x; s) and G-ensemble for the other). 

A dissymmetry between the two is introduced by formulating a necessary 
and sufficient condition for the stability of equilibrium configurations of one 
ensemble, say the F-ensemble: 

d 2 x F(x = X l (s);s)<0 (5) 

(i. e. that F be a maximum). Thus any change of stability, in the F-ensemble, 
will occur if and only if d^F crosses zero. 

^or the reader who wants to have in mind a specific example, s may be interpreted 
as the total energy E (then by working at fixed s one has the microcanonical ensemble), 
the variable x viewed as the energy repartition between two subsystems and the potential 
F(x; s) may be interpreted as the entropy out of equilibrium when the variable x is not 
at its equilibrium value Xi(s). Then t is the inverse temperature and the potential — G is 
the free energy divided by the temperature. The reader might usefully consult ref. 19 to 
see in detail how this method apply to situations containing a black hole surrounded by 
radiation. 
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Contrariwise, in the G-ensemble, at fixed t, s fluctuates^] as well around 
the solution of eq(2): S(x,t). Thus further conditions must be met to fulfill 
stability (i.e. G has to be a true maximum and not only a saddle point 
in order to insure that the new fluctuations be bounded). Since those new 
conditions need not be met even when eq(5) is satisfied, the G-ensemble is 
the less stable one. Our purpose is now to prove that those new conditions 
inevitably restrict the range of stability of equilibrium situations in the G- 
ensemble. (At the end of the paper, we shall discuss the "opposite" problem, 
that is, how, when one's consideration begins with the knowledge of the 
stability limits of the less stable ensemble, one may compute the extension 
of the stability range for the more stable ensemble.) 

We will first review how changes of stability, in the F-ensemble, manifest 
themselves using Poincare method. The reader unfamiliar with this method 
may usefully consult refs [1,2,19]. In addition, we are using the same no- 
tations and conventions. The reason for this review is that it affords an 
instructive illustrative example of how the topological behavior of the linear 
series in the vicinity of a critical point inevitably restricts the stability range 
of the Legendre transformed ensemble. 

The idea of the method is to construct the linear series, locus of solutions 
of eq(3), as a function of s : x = Xi(s). Poincare theorem states that stability 
changes may occur only at a bifurcation (i.e. when two series intersect) or 
at a turning point (i.e. when two series merge into each other). Bifurcations 
occur because of an "excess" of symmetry and can be removed by the slightest 
modification of the potential F(x; s)t 2 ' 16 '. Hence, only the behavior of the 
series near a turning point is studied here. 

More specifically we consider the particular situation where one has two 
equilibrium solutions of eq(3) for s smaller than a certain maximum denoted 
Sf. These solutions merge into each other as s tends to s/ and do not exist 
for s greater than Sf. (This situation corresponds to the case i (Fig. 2) of 

2 We are using the thermodynamical terminology even though the present analysis ap- 
plies as well to mechanical stability, thus by 'fluctuates' we mean that we have to consider 
neighboring solutions around the equilibrium one in ordre to determine stability. In the 
same spirit by 'fluctuations of x 1 we shall designate the quantity (— d^F)^ 1 / 2 which is the 
RMS fluctuations of ir in a statistical ensemble weighted by e F . We refer again to ref. 19 
for explicit examples in which the notions of stability and fluctuations are defined more 
concreately. 
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the classification of ref.l) The series of solutions x = X\(s) describes stable 
equilibrium situations for which eq(5) is satisfied, and the series x = X 2 (s) 
describes unstable situations. The method then proceeds by evaluating the 
second total derivative along those series: 

^ = <?F . +aAF .«. (6) 

where F a = F(X a (s), s) and a = 1,2. With the help of the following identity 

^ (d x F a ) = = d s d x F a + d 2 x F a ^ (7) 

eq(6) becomes 

^ = dlF a - (d x d s F a ) 2 /dlF a (8) 

Hence the change of stability at s — Sf, where d x F a = 0, is manifested by 
the divergence of this second derivative unless 

d x d s F a (s = s f ) = (9) 

(We do not consider this degenerate possibility, which leads to a bifurcation, 
since it has been provent 2 ' 16 ! that the slightest modification of F will destroy 
the simultaneous vanishing of d 2 F a and d s d x F a ). We emphasize the role and 
the physical meaning of eq(7). If d x d s F a ^ 0, one has 



dX a 



{dlFaY 1 (10) 



ds 

which means that near the instability, the evolution of the equilibrium value 
X a (s) (or the mean value in thermodynamical ensembles) is entirely con- 
trolled by the growing of the fluctuations. Hence, the knowledge of X a (s) 
is sufficient to determine the stability changes without computing separately 
the fluctuations themselves. 
Furthermore since 

d 2F a (^(tfFY 1 ^ +oo along X x (s) 

ds 2 s^f { ) \°* ra ) _> _oo along X 2 (s) 1 J 

the linear series 1, describing stable situations, lies inevitably below (in the 
s — t plane) the series 2, as depicted in Fig.l. 
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Before examining what are the consequences of this unique behavior, we 
briefly discuss what could have happened in more complex situations. 



1. For ensembles with many x variables, the linear series 1 could describe 
already unstable situations with n negative modes (where n < number 
of degrees of freedom described by the x variables). Nevertheless by 
going anti clock wise along the curve of Fig.l. one still increases the 
number of unstable modes each time d s F a has a vertical tangent. 

2. In the cases where one has many parameters s, each of them or any non 
singular function of them can be used to study the changes of stability 
of the F-ensemble, leading to identical predictions because the vertical 
tangent only occurs when d x F vanishes. But if one wants to study the 
stability of a particular Legendre transformed ensemble and compare 
it to the stability of the initial ensemble, it is appropriate to select the 
parameter which is used to define the Legendre transformation. (In the 
present case with only one parameter, we are automatically in such a 
situation.) 

We now examine the stability in the G-ensemble. We first redefine in a 
slightly different manner the Legendre transformation, eq(l), to see clearly 
how the transformation from F to G enlarges the spectrum of fluctuations 
by one eigenvalue^ (in our case we have two fluctuating quantities x and s) 
and why the new fluctuations always encode the changes of stability of the 
G-ensemble. 

G(x,s;t) = F(x;s) - st (12) 
Equilibrium configurations are now provided by 



d x G = = d x F (13) 
d s G = = d s F-t (14) 

Eq(13) leads identically to eq(3) since now the only x dependence of G is 
through F. Eq(14) connects s with its Legendre-conjugate t as in eq(2). 

3 We recall that in a canonical ensemble the total energy E (here s) of the system 
fluctuates due to the contact with the heat-reservoir. 
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These two equations furnish: x = X a (t) = X a (S a {t)), where S a (t) are the 
two values of the solution of eq(14) evaluated at the equilibrium situations 
x = X a . 

The range of stability is obtained as previously by computing the second 
derivative of G, with respect to t, along the linear series x = X a (t), see 
eqs(6-ll) 

^-ZG. + W.^ + BAG.f (15) 
where G a = G(X a (t), S a (t);t). Since by construction in eq(12) one has 

d?G = 0, d t d s G = 0, d s d t G = -l (16) 

one finds 



d 2 G a dS a 
dt 2 ~ dt 

= -(dlF a -{d x d s F a f/d 2 x F a y l (17) 

(d 2 F a y l _ (dTA- 1 
\ds 2 ) ~ ^ ds ) 

where T a (s) = d s F a and where the second equality follows, as for the F- 
ensemble (see eq(7)), from 



d (d x G a ) = = %fJ%- + d s d x F a dSl ' 



dt dt dt 

4 (d s G a ) = = -1 + d x d s F a ^ + d 2 F a ^ (18) 
at at at 

We thereby recover the well-known identity, eq(17), between first total deriva- 
tives. The reasons that we prove this identity are the following: 

1. By the very structure of the Legendre transformation which implies 
eqs(16), we see that the fluctuations of s ( i.e. (—d 2 F/ds 2 )~ 1 / 2 ) alone 
control the vertical tangents of d t G a = S a (t). But they also control, as 
we will prove, the changes of stability in the G-ensemble. 
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2. These new fluctuations encode in a very specific way the previous fluc- 
tuations of x (see eq(17)) and one has still to impose eq(5) in order 
to have stable equilibrium (even though, there will be no new unsta- 
ble mode, in the G-ensemble, when d 2 F a will become negative, see the 
Appendix). 

From eq(17) we obtain that changes of stability in the G-ensemble may 
occur when d 2 F a /d 2 s = and that (since only maxima describe stable equi- 
librium) the more stable series is the one with d 2 F a /d 2 s < 0, i.e. see eq(8) 
when: 

d 2 F a <(d x d s F a ) 2 /d 2 x F a (19) 

Since d 2 F a is negative along the stable linear series 1 and tends to — oo as s 
tends to s/, the inequality (19) has to break down before sj at some maximum 
point s = s g . (We assume that F has regular second derivatives). Hence, 
since d 2 F a /ds 2 crosses zero before d 2 F a) the new fluctuations will always 
control the changes of stability of G. This proves the third and the last point 
mentioned in the introduction. 

One can visualize this ordering by examining the behavior of the linear 
series in the vicinity of sj and s g . Only two situations may occur. Either 
eq(19) is never satisfied and the G-ensemble is always unstable with one 
negative mode (case I of Fig.2. An example of this case is provided by a 
Black Hole in contact with a heat bathJ 7 !), or eq(19) is satisfied until s = s g 
and the G-ensemble is stable up to that point (case II, see the example of 
star cluster in ref. 1). 

For the interested reader we present in the appendix the (more tradi- 
tional) analysis of the stability of the G-ensemble when any reference to s 
has been eliminated, through the use of eq(2), before looking at equilibrium. 
The result one obtains is that the quadratic fluctuations of x now evaluated 
at fixed t, rather than fixed s, encode automatically the determinant factor 
d 2 F a /ds 2 hence leading again to the same range of stability. In addition 
we present in this appendix how, starting from G, one recovers the correct 
enlarged range of stability for the F-ensemble by an inverse Legendre trans- 
formation. The mathematical reason for which one increases the stability 
range, by going back from G to F, is directly related to the fact that when 
both ensemble are stable one has 
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d 2 F c 



d 2 G, 



d 2 F c 



) 



i 



a 



<0 , 



a 



a 



> 



(20) 



ds 2 



dt 2 



ds 2 



This inevitable discrepancy of signs has the following physical interpreta- 
tion. Parameters of the s type, like energy, angular momentum or volume 
represent conserved quantities of the Hamiltonian of an isolated system and 
have to be treated as the variables x (stability upon extremization, eq(14), 
requires the same sign for the curvature, see eq(19)). In contrast parame- 
ters of the t type, like temperature, angular velocity or pressure, are kept 
fixed by the intervention of an external "reservoir" . Furthermore they enter 
only in a linear way into the state function (see eq(12)) and play an identi- 
cal role to Lagrange multipliers. The fact that their extrema correspond to 
minima, eq(20), translates in Statistical Mechanics into the fact that inverse 
Laplace transformations behave like Fourier transformations' 15 '. Thus their 
integration, due to their linear dependence, leads to a Dirac 5-function which 
reduces by one the fluctuating quantities. 

In thermodynamics, in complex situations, one deals with many control 
parameters. The ensemble with all its parameters being of the s type is the 
most stable one and describes the most isolated system. The stability of this 
ensemble does not refer to any other ensemble and the extremization of its 
state function is a rephrasing of the Second Principle of Thermodynamics^. 
Contrariwise the stability of all other ensembles is subject not only to addi- 
tional necessary conditions (like the positivity of the specific heat, see eq(19)) 
but also to the stability of the most stable ensemble, eq(5). 

Application to rotating black holes 
To illustrate the power and the simplicity of our analysis we consider the 
situation in which one has a rotating black hole in contact with a heat bath. 
Working at fixed angular momentum, it has been noted that there exists 
a temperature at which the specific heat changes sign' 8 ' 9,10 ]. The question 
is then: is the equilibrium, in the canonical ensemble, stable or only less 
unstable when the specific heat is positive? 

The problem of the stability of a rotating black hole was addressed in 
ref(10). In the microcanonical ensemble (i.e. at fixed mass M, fixed angular 
momentum J and with the entropy S as the relevant state function), it was 

4 This was pointed out to me by R. Brout. 
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stated that the absence of vertical tangent in the plot (3 = OmS versus M 
indicates that their is no change of stability when one varies the ratio J/M 2 
(see Fig. 3). Then, since it had been proven that isolated Schwarzschild 
black holes are stable [17] , it was correctly concluded that all isolated rotating 
black holes are stable^. 

In the canonical ensemble, by virtue of this analysis and our statement 
2 (see Introduction), one immediately concludes that the flip of sign of the 
specific heat at (3 = (3 C (see Fig. 3) does correspond to a real change of 
stability and that all equilibrium situations along the linear series with the 
smaller mass are stable. This conclusion is stronger that the one given in 
ref.(10) because the authors did not take into account the fact that the 
instability of the Schwarzschild black hole in contact with a heat bath is due 
to a single instable mode. 

Appendix A 

Two points will be discussed in this appendix. The first one is the recover- 
ing of the changes of stability of the G-ensemble when all dependence of s 
has been eliminated prior looking at equilibrium configurations. The second 
point concerns the extension of the stability range when, starting from G, 
one wishes to analyze the stability of the F-ensemble. 

Since all explicit dependence of s has been eliminated, one has to ana- 
lyze the stability of x fluctuations at fixed t rather than fixed s. We have 
nevertheless to verify that the fluctuations of s at fixed x described by <9 2 F 
are bound {i.e. d^F < 0. We assume here that this is satisfied). The state 
function which control the x fluctuations is (see eq(l)): 

G(x; t) = F(x; S(s; t)) - tS(x; t) (A. 1) 

where S(x;t) is the solution of eq(2). Equilibrium situations are furnished 

(see eqs(3,4)) by: 

Nevertheless, since the entropy S(M, J) is already an extremized quantity with respect 
to the unknown variables "x" which describe the microscopical states of the black hole, 
one has to assume that a simultaneous vanishing (see eq(9) and the associated discussion) 
which could hide changes of stability does not occur. 
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d x G\ t = = d x F\ a (A. 2) 

by virtue of eq(2). (In this appendix, in order to prevent any quibble we 
specify by a vertical line which quantity is kept fixed upon derivation.) 
Stable equilibrium requires G to be a maximum (see eqs(5, 17,19)) 

dlG a \ t < (A. 3) 

This condition reads 

d 2 x G a \ t = d 2 x F a \ s + d s (d x F a \ s )\ x (dSJdx) \ t (A. 4) 

where S a (t) = S(X a (t),t). (dS a /dx) \ t ) can be obtained by taking the total 
differential of eq(2): 

d x (d s F\ x )\ s dx + d 2 s F\ x ds = dt = (A. 5) 



hence 

dx 

Then eq(A.4) becomes 



<1S \ = -(d x d s F)/(d 2 s F) 



d 2 x G a = d 2 x F a -(d x d s F a ) 2 /d 2 F a 

= ( A - 6 ) 

where we have used eq(8). By virtue of the analysis which follows eq(19) and 
which indicates that d 2 F a /ds 2 crosses always zero before d^F a , one recovers 
that the changes of stability of the G-ensemble occur when d 2 F a /ds 2 van- 
ishes. Furthermore, eq(A.6) indicates that the vanishing of d 2 F a at s — Sf 
does not lead to the vanishing of d 2 G a because of the divergence of d 2 F a /ds 2 
at that point (see eq(ll)). Finally, by an analysis similar to the one given 
after eq(18) one can easily show that the second negative eigenvalue of the 
G-ensemble will always appear after the first one of the F-ensemble, thereby 
reinforcing the point 1 (see Introduction) by introducing a well defined or- 
dering. 
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The second point of this appendix concerns the following problem: having 
obtained the limit of the stable configurations for the less stable ensemble 
(when A. 6 vanishes) how can one evaluate the extended range of stable config- 
urations in the F-ensemble? This problem is by no means purely academical: 
Under certain circumstances the evaluation of quadratic fluctuations can be 
performed easily only in, say, the grand canonical ensemble. Then two ques- 
tions arise when one wants to use those results into the evaluation of the 
stability limits for more stable ensembles (see refs. [4,5,15], where this prob- 
lem was encountered and discussed). The first one concerns the possibility 
of extending the evaluation of the fluctuations outside the stability domain 
of the ensemble in which they were computed. The second question is: how 
the instability of the initial ensemble will manifest itself into the expression 
of the stability of the fluctuations in a more stable ensemble? Answering 
those questions can be achieved in many different ways. We only sketch two 
possible ways using our simple mathematical abstraction. 

The most pedestrian way takes the lines of the analysis presented from 
eq(12) to eq(19): one introduces a third fluctuating quantity (i.e. t) and one 
finds that the role of its fluctuations is to suppress completely the fluctua- 
tions of s thereby leading back to the "initial" F-ensemble with its single 
fluctuating variable: x. (This procedure corresponds to the second (noted 
II) diagonalization scheme of ref. [15]). 

A second illuminating method proceeds, as shown above (see eq(A.6)), 
by expressing the quadratic fluctuations of x at fixed s in terms of the fluc- 
tuations evaluated at fixed t. One has: 



Suppose now that one increases s, starting from stable configurations for 
G (i.e. for s < s g , see Fig. 2, where d 2 F a < 0). When s crosses s g both 
d 2 G a and d 2 G a /dt 2 flip sign but their product is perfectly well defined. Hence 
d 2 F a remains negative and will vanish only when d 2 G a /dt 2 vanishes. Eq(A.7) 
proves that the evaluation of the x fluctuations in the F-ensemble (controlled 
by d 2 F a ) can be computed safely in the G-ensemble even in situations which 
correspond to unstable equilibrium in G. (This simultaneous vanishing, at 
s = s g , is precisely what was encountered in ref. 4, and subsequently analyzed 
in refs. 5 and 15. It leads to the first (I) diagonalization scheme of ref. 15.) 




(A. 7) 
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Figure captions 



Fig. 1. Plot of the equilibrium curve T a (s) = d s F a (s), evaluated along the 
linear series Xi(s) and -^(s), in the vicinity of the turning point s = s g . 
The series Xi (s) describe stable equilibrium configurations and lie inevitably 
under the series X 2 (s). At s = s g , the fluctuations of x are unbound {d^F a 
vanishes) and this manifests itself by a vertical tangent of T a (s). This can be 
understood by extending the definition of the parameter t away from equi- 
librium: t(x, s) = d s F, and letting it fluctuate around T a (s) according to the 
x fluctuations' 19 ^. 

Fig. 2. This is the same plot as the one of Fig. 1. The two possible behaviors 
of T a (s) are displayed. Along the linear series (1,1) the G-ensemble is always 
unstable with one negative eigenvalue. Along (1,11) the G-ensemble is stable 
for s < s g (or t < t g ). The horizontal tangent of dS a /dt at t = t g , indicates 
that the fluctuations of the variable s are infinite. This should be contrasted 
with the parametrical fluctuations of t in the F-ensemble. 

Fig. 3. Plot of the equilibrium curve (3 = 8mS(J, M) versus the mass M, at 
fixed angular momentum J and describing a Kerr black hole' 7 ' 9 l At (3 — (3 C , 
8mP flips sign. Below M c , in the canonical ensemble, the black hole is in 
stable equilibrium with the heat bath. This situation corresponds to the 
case II of Fig. 2, with s/ rejected at infinity. 
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